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Abstract 

This paper gives a review of our recent work on the quantum Frenkel-Kontorova 
model. Using the squeezed state theory and the quantum Monte Carlo method, we 
have studied the effects of quantum fluctuations on the Aubry transition and the 
behavior of the ground state wave function. We found that quantum fluctuations 
renormalize the sinusoidal standard map to a sawtooth map. Although quantum 
fluctuations have smeared the Aubry transition, the remnants of this transition are 
still discernible. The ground state wave function also changes from an extended 
state to a localized state. The squeezed state results agree very well with those from 
the Monte Carlo and mean field studies. 

Key words: Frenkel-Kontorova model; squeezed state; quantum phase transitions; 
quantum chaos 
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1 Introduction 

Many physical systems sucli as cliarge density waves, magnetic spirals, Joseph- 
son junctions and adsorbed monolayers exhibit a competition between two or 
more length scales. To understand such modulated structures, a simple model, 
known as the Frenkel-Kontorova (FK) model [1], was introduced more than 
half a century ago. The FK model describes a chain of atoms connected by 
springs in the presence of an external sinusoidal potential. It has two length 
scales: the natural length of the spring and the period of the external potential. 
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Due to the competition between these two length scales, the FK model ex- 
hibits a wealth of interesting and complex phenomena. Extensive studies of the 
FK model have been made since its introduction. In most of the earlier stud- 
ies, the FK model was treated in the continuum approximation. Although the 
continuum approximation, which reproduces the sine-Gordon equation and its 
soliton solutions, gives some very useful information, it misses many essential 
features. A major advance in the study of the FK model was made by Aubry 
[2]. Aubry reverted to the discrete version of the FK model and made essential 
use of the seminal Kolmogorov-Arnol'd-Moser (KAM) theorem. In analogy to 
the breakup of a KAM torus, Aubry discovered the "transition by breaking 
of analyticity." The Aubry transition is similar to a phase transition, thus all 
the apparatus invented in the study of second-order phase transitions can be 
borrowed to study this transition. A fairly complete picture of the FK model 
was thus obtained. 

So far most of the studies made on the FK model have been confined to the 
classical model. Very little study has been done on the quantum FK model. 
Although the FK model looks simple, it's nevertheless a non-integrable inter- 
acting many-body system. There aren't that many methods at our disposal 
to study the quantum behavior of such a model. Moreover, besides theoretical 
interest, there are also practical considerations that demand a deeper under- 
standing of the quantum FK model. For example, the study of friction, or 
tribology [3], is a very important subject both in terms of fundamental sci- 
ence and applications. And indeed, the FK model has been used to study 
tribology. However, as technology is moving into the nano-regime, nanotribol- 
ogy [4] is also becoming a very important subject of study. In the nano-regime, 
it's more likely than not that quantum mechanics will play an important role. 
This makes an understanding of the quantum FK model not only desirable but 
very likely imperative. This is the motivation for our interest in the quantum 
FK model. 

The first study of the quantum FK model was made by Borgonovi, Guarneri 
and Shepelyansky [5]. They mainly used a numerical approach. Later, Berman, 
Bulgakov and Campbell [6] studied the same problem using a mean-field ap- 
proach. In our study we will use a hybrid analytical and numerical approach. 
The analytical method is based on the squeezed state theory, and the numer- 
ical method is based on the quantum Monte Carlo (QMC) algorithm. 

This paper gives a review of our recent work on the quantum FK model [7,8] . It 
is organized as follows. In Sec. 2, we will define the quantum Frenkel-Kontorova 

model and use the Monte-Carlo method to study it. In Sec. 3 we will formulate 
the squeezed state approach. In Sec. 4 we will study the quantum effects on 
the Aubry transition. In Sec. 5 some concluding remarks will be given. 
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2 Quantum Frenkel-Kontorova model 



The Hamiltonian of the quantum Prenkel-Kontorova model is described by: 



^2 

^ + -{xi+i -Xi- of - y cos(go^i) 



(1) 



Here m is the mass of the particles, 7 the elastic constant of the spring, 271 /q^ 
the period of the external potential, V the magnitude of the external potential, 
and a the equilibrium distance between two nearest neighbor particles. For 
convenience, we shall rescale the variables so as to obtain a dimensionless 
Hamiltonian 



P2 1 . 



Xi — nY — K COS Xi 



(2) 



where Xi = qoXi, Pi 



"effective Planck's constant" : 



= qoa, and K = Vq^/j. We define an 



h 



n. 



(3) 



This effective Planck's constant is the ratio of the natural quantum energy 
scale TiUjQ to the natural classical energy scale 7/^0 > ^^id — 7/m. This 
quantity provides a measure of the importance of quantum effects. In par- 
ticular, quantum effects become crucial at zero temperature when thermal 
fluctuations vanish. 



In the classical FK model, the positions of the particles in the ground state, 
denoted by X?, are determined by the following equation: 



Xt,,-2Xt + XU = Km.Xt. 



(4) 



This equation can be cast into the form of the standard map by defining 



h+i^ li + K Sin Xf, 



i+1 + ^i- 



(5) 



To describe the transition from the unpinned phase in the subcritical regime 
to the pinned phase in the supercritical regime, the hull function is introduced. 
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It is defined by 



(6) 



Here Xi is the unperturbed {K = 0) particle position. Another useful function 
is called the gf-function, defined by 

9 = ^{Xt^i - + Xl^ (7) 

For K < Kc, the ground state of the classical FK model corresponds to an 
invariant curve in the standard map (5). The hull function is a monotonic 
analytical function [Fig. 1(a)], and the gi-function is g{X) = sinX [Fig. 1(b)]. 
At K — Kc, the last invariant curve in the standard map breaks up and 
becomes a Cantorus. The hull functions becomes a monotonic function with 
a countable number of discontinuities [Fig. 1(c)], and the g^-function becomes 
a Cantor subset of the function sinX [Fig. 1(d)]. 

By using the quantum Monte Carlo method, we can calculate the ground state 
averages of the particle positions Xi and plot the following three functions so as 
to compare with their classical counterparts: (1) The quantum hull function, 
i.e., Xi as a function of the unperturbed positions Xi, (2) the quantum g- 
function which is defined as 

g,^^(Xi+^-2Xi + Xi_,); (8) 

and (3) {Xi, li) in the phase space of the standard map. 

The quantum effect does not destroy the functional dependence of g^ on Xi, 
instead it renormalizes it to a sawtooth map [Fig. 2]. This phenomenon was 
first observed by Borgonovi et al. [5] and was interpreted as a consequence of 
quantum tunneling between the edges of the gap in the hull function. 



3 Squeezed state approach 

In this section, we will study the effect of quantum fluctuations on the FK 
model by using the squeezed state approach [7]. 

The squeezed state was first used in quantum optics and has proven to be very 
useful in many other fields [9-13]. The squeezed state approach starts from 
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the time-dependent variational principle 

S I dt{^{t)\th--H\^{t))^0. (9) 

Variation with respect to (w.r.t.) gives rise to the Schrodingcr equation. 

The exact solution may be approximated by restricting the choice of states 
to a subspacc of the full Hilbcrt space. In the squeezed state approach, the 
squeezed state is chosen as |^l'(t)). 

The squeezed state |$(q;, is defined by 

= e^(")e^(^)|0), (10) 

where 

-^(q^) = - (11) 



is the Weyl's displacement operator, and 

1 



f{f3) = -j:{a^Aja]^d,Pldj). (12) 



is the squeezed state operator. S^{a) = —S{a),f\(3) = — f{P). For brevitiy, 
we will use the abbreviation |$) = |$(a,/9)). If = 0, the squeezed state 
is reduced to the coherent state. However, as we shall see later, the coherent 
state is not adequate to study quantum fluctuations. 

The position and momentum operators for the ith particle are written as 



U t ^ ' 

Pi = iVf(a- -di). 



Here, d\ and di are the creation and annihilation operators which satisfy the 
canonical commutation relations: [Sj, a]] = ^jj, [cj, dj] — and [a|, a]] = 0. 

Using 1$) as a trial wave function, we can flnd the expectation values of the 
position and momentum operators of the ith particle [9] : 

X, = ($|X,|$) = y|« + a,), ^^^^ 
P, = {^P^\^)^-i^{a*-ai). 
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Fluctuations in the position and the momentum are given by 



AX2 = ($|(X,-X,)'|$) = ^Gn, 

AP^ = ($1 (P, - P,)2|$) = h + 4 J2 ^uGikIlk}j . (15) 

The fluctuation covariance between the ith particle and the jth particle is 

AXiAXj = - Xi){Xj - Xj)\^) = hGij, (16) 



where 



G,, = Kcosh^ ^f3f3t + sinh^ ^(3(3% + \{M(3 + (3^M)ij, ^^^^ 
n,,- = iGr.i(M/3-/3tM),,-, 



and, 



^ sinh^/^cosh^^ 



Since /3 is a symmetric matrix, G^- = Gji and Iljj = Iljj. Using the following 
relation 

($1 cosX,|$) = exp (-^G.A cosXi, (19) 



we obtain the expectation value of the Hamiltonian, H = {^\H\^), 



i y l,k j i 

+ E ^ + ^mi+i - SG.+h) - X^i^exp [~G^ cosXi. (20) 



The variables Xi and Pj, and Gj^ and 11 form explicitly canonical conju- 
gates [10]. To find the ground state of the quantum FK model, we will use a 
variational approach in which these four variables are treated as independent 
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variables. Variation with respect to Pi and Xi yields, respectively, 

Pi = 0, 

_ _ _ _ (21) 

Xi+i - 2Xi + Xi_i = Ki sin Xi, 



where 



Ki^Kexp(-^G^A. (22) 



This equation determines the expectation value of the particle's position. Un- 
like its classical counterpart {h = 0, Ki = K), this equation is coupled to the 
quantum fluctuation by hOu. Variation w.r.t. Uij leads to 

^GikUkj^O. (23) 

k 



To obtain the equation for Gij, we first take variation w.r.t. Gik and note that 

^ = (24) 



We then multiply both sides of the equation by Gkj and sum over k. We then 
get the closed set of equations for the covariance Gij: 

{GF)ij = Gi-ij + Gi+ij, (25) 



where 

F., = (l + f cosX.)^,,-^^. (26) 



This is a set of equations determining the quantum fluctuations of the par- 
ticles. G is an X symmetric matrix which provides information about 
fluctuations. Its diagonal elements give the variance of each particle while its 
off-diagonal elements give the covariance between particles. Prom this infor- 
mation we can calculate the correlation function of the quantum fluctuations. 
These equations are coupled to the expectation values Xi. 

Up to this point, we have obtained |iV x (A^ -|- 1) -|- equations for all the 
variables. These equations provide some qualitative information about the 
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system. If we define li+i — Xj+i — Xi, Eq. (21) can be cast into tfie form of a 
standard map, 

Jj+i = Ii + Ki sin Xi, 

- - - (2' J 

Xi+i — /j+i + Xi. 



Similarly, if we define Qi+ij — Gi+ij — Gij, Eq.(25) can also be cast into the 
form of a map: 

Qi+ij = Qij + [G{F - 2)]ij, 
Gi+ij = Gij + Qi+ij- 



The difference between the classical case {h = 0) and the quantum case is 
readily seen in Eq.(27) . In the classical case the coupling constant does not 
change from potential to potential; however, in the quantum case, the effective 
coupling constant changes from particle to particle. Because Gu > for any 
h ^ 0, Ki < K. which implies that the quantum fluctuation will reduce the 
coupling strength. Another important difference is that in the classical case the 
positions of the particles in the ground state are determined by the standard 
map whereas in the quantum case they are determined by (iV + l)-coupled 
two-dimensional maps. This makes the quantum FK model extremely difficult 
to study analytically. 

One notices that when /5 = 0, Gij = 1/2, which is the result of the coherent 
state theory. This is not the case in the quantum FK model. In Fig. 4 we plot 
hGij as a function of i and j. The data are obtained from QMC. It is obvious 
that Gij 7^ 1/2. Thus it seems that the coherent state is not adequate for the 
study of the quantum FK model. 

We now make some comparisons with the QMC results. As mentioned before, 
to find the solution of these two sets of equations [Eqs.(21) and (25)] is equiva- 
lent to finding a periodic orbit in the 2(A^ + l)-dimensional map [Eqs.(27) and 
(28)]. This is a very difficult problem to solve even numerically. Nevertheless, 
we will try to solve Eq.(21) self-consistently to see whether this equation can 
indeed give rise to the sawtooth map [5] . 

Using QMC, we obtained the expectation value of the particle's position Xj, 
from which we can construct the quantum hull function /g, as shown in Fig. 
2(a) for K — 0.5 and (c) for K — b, respectively. The quantum g^g-function is 
shown in Fig. 2(b) and (d) for K = 0.5 and K = 5, respectively. The quantum 
fluctuation hGn calculated from QMC is shown in Fig. 3 for the case K — 5. 

As usual, we use Q particles and P external potentials with period 27r, and the 
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periodic boundary condition is imposed [14]: Xq+i — Xi + 27rP. The results 
shown in the figure are for P/Q — 34/55. 

To compare the squeezed state results with those from QMC, we substitute Gu 
calculated from QMC into Eq.(21). We then compute the expectation value of 
the particles' positions by using Aubry's gradient method. We then construct 
the quantum hull function fq and the quantum (yfg-function shown in Fig. 5. 
The results agree surprisingly well with those from QMC. The most striking 
feature is the sawtooth shape of the g^q-function. This phenomenon was first 
observed by Borgonovi et al. [5] in their QMC computation and has been 
explained as a tunneling effect. Later, Berman et al. [6] obtained this result 
by using a mean field theory. In the framework of the squeezed state theory, 
the renormalization of the standard map to a sawtooth map is a consequence 
of quantum fluctuations. Our results show that the squeezed state approach 
indeed captures the effects of quantum fluctuations. 



4 Quantum effects on the Aubry transition 

In this section we will study the effects of quantum fluctuations on the Aubry 
transition. We will flrst study the behavior of the wave function of an incom- 
mensurate ground state. We will then introduce various indicators that might 
be used to characterize the quantum Aubry transition. 

4..1 Ground state wave function of the quantum Frenkel-Kontorova model 

We will use the QMC [15] to study the ground state of Hamiltonian (2). As 
in the classical case, we will concentrate our attention on the incommensurate 
state characterized by the golden mean winding number ujg — ("n/S — l)/2. 
In the classical case the winding number is defined as the average separation 
of atoms per period, i.e., u = limjv^oo ^gAf^f " • H is a rational number, the 
ground state is commensurate; and if it is an irrational number, the ground 
state is incommensurate. In the quantum case, Xi is defined as the expectation 
value of the position of the ith particle. As usual, we use the method of 
continued fraction expansion and approximate Ug by its rational convergents 
i^n — Fn/Fn+i, whcrc F„ is the nth Fibonacci number defined by the recursion 
relation F^+i = + with Fq = 0, Fi = 1. In our computation, we choose 
Fn+i particles embedded in F„ external potentials. 

Since the external potential is periodic with period 2tt, we can fold the wave 
function to this period and then take the average over all the particles in the 
interval [0,27r]. This quantity gives the probability of finding the particle at 
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a given position X. Wc plot the averaged probability density in Fig. 6 

with 144 particles embedded in 89 potentials for a fixed h (= 0.2) but different 
values oi K. Here is normalized, i.e., J^'" dX{\'^{X)\'^) = 1. We observe 

that, in the small K regime, the probability of finding the particles at any point 
of the potential is almost the same (see, for example, the curve for K — 0.1). 
This is quite similar to the unpinned phase in the classical case. However, as 
the coupling constant increases, the probability of finding the particle in the 
upper part of the potential decreases whereas the probability of finding the 
particle in the lower part of the potential increases, as is noticeable from the 
appearance of peaks in the curves. As the coupling constant increases, the 
probability of finding the particles at the top is almost zero (see, for example, 
the curves for = 2 and 5). 



4-2 Indicators of the Auhry transition 



4-2.1 Disorder Parameter 

In the classical FK model, Coppersmith and Fisher [17] have proposed a "dis- 
order parameter" to describe the transition from the pinned phase to the 
unpinned phase. This parameter is defined as the minimum distance of a par- 
ticle from the top of the potential, D — min^^^ \Xf — 27r(n + |)|. If 7^ 0, 
the particles are pinned; when D = 0, they are unpinned. This "disorder pa- 
rameter" also measures the discontinuity (or width of the biggest gap) of the 
hull function. One might want to use the same function to describe the quan- 
tum crossover. For instance, one may define a very similar quantum disorder 
parameter Dg — mini_„ |Xj — 27r(n -|- |)|, where Xi is the expectation value 
of the ith particle's position. However, this parameter Dg could not capture 
the crossover of the ground state wave function. The reason is that, in the 
quantum case, the particle can tunnel from one side of the potential to the 
other. Thus the gap in the classical hull function does not survive quantum 
fiuctuations (see Refs. [5] and [7]). It turns out that for large K, Dg might 
still be close to (or equal to) zero. Therefore, a new parameter is needed to 
describe the wave function crossover. To this end, we define the probability of 
finding the particles at the potential top as a quantum "disorder parameter" , 



N-l ^ r I 



X-27r(n + -) 



dX. (29) 



In Figs. 7(a)- (b), we plot Pt as a function of K for different temperatures Tg 
(Te = h/r, where r is the length of the imaginary time axis) and different 
system sizes. In Fig. 7(a), we fix h{= 0.2), ut = 13/21 and set the temperature 
Te = 0.2/30,0.2/120,, 0.2/480. The convergence is quite fast as is seen from 
the figure. The sharp decrease of Pt is very evident for 1 < X < 2. In this small 
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K regime, Pt changes significantly: it drops almost two orders of magnitude. 
This drastic change can also be seen from Fig. 7(b). There the temperature is 
fixed at = 0.2/120, which can be regarded as effectively zero temperature. 
The system size is changed from 21 particles to 144 particles. 



^.2.2 Maximum Variance 

Another parameter that can be used to depict this crossover is the maximal 
fluctuation of the particles. In the quantum case, we have observed that the 
particle situated at a position nearest to the top of the potential always has 
a maximal fluctuation since it has a higher classical potential energy. This 
observation has been verified numerically by QMC and theoretically by the 
squeezed state theory [7] . Thus we can use this maximal fluctuation as another 
measure of quantum "disorder" . To make a qualitative comparison with the 
classical disorder parameter we define 



Here the average (■ ■ ■) is taken over all the paths (approximately 4,000) used 
in QMC. A as a function of K is plotted in Fig. 8. The computations given 
in this figure have been carried out with uj^ — = 34/55. 

It is noticeable that the transition of the wave function is characterized by the 
different fT-dependent behavior of A. In the small K regime, A is a constant 
independent of K but changes with h. In the large K regime, A increases 
with K but does not change with fi. Furthermore, A is approximately equal 
to the classical disorder parameter D. For comparison, we have included D in 
the inset of Figure 8. In the classical case, the transition of course occurs at 



These results show that the width and the shape of the probability density 
do not change with K in the small K regime. The width only depends on the 
strength of the quantum fluctuation Ti. However, this picture changes signifi- 
cantly for large K. In this regime, the profile of the probability density spreads 
out, and the width of the probability density is quite insensitive to quantum 
fluctuations and depends only on the coupling constant. It should however be 
stressed that the analogy between D and A cannot be carried too far since A 
is not exactly the quantum correspondence of D. This is why even if we let fi 
go to zero, A would not be zero in the small K regime. 




(30) 
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4-2.3 Correlation function 

As another indication of the quantum signature of the Aubry transition, we 
study the correlation function between the particles' fluctuations. The cross- 
correlation function Cy is defined by 



It describes the fluctuation correlation between the particles i and j. lij = 
{{Xi — Xi){Xj — Xj)) is the covariance of the fluctuations between the ith 
particle and the jth particle. Here (■ • •) denotes the ensemble average over 
space and 'time' in QMC [14]. Since lij is symmetric, Cij — Cji. In our nu- 
merical calculations, we take 4000 paths and each path is divided into many 
steps in time, ranging from 1500 to 3000 depending on h. The convergence is 
guaranteed by increasing the number of steps and the number of paths until 
the same results are reached. 

Figure 9 shows the contour plots of Crj for K = 0.5 and K = b. The plots 
are drawn in 12 equal steps from —0.15 to 1. A negative correlation function 
implies that the particles' quantum fluctuations are out of phase. The sharp 
change of Cij from K — Q.b to K — b \s very readily noted. In the large K 
regime, the correlation function matrix C^j is almost diagonal with a very small 
band width whereas in the small K regime the number of nonzero off-diagonal 
elements increases and the band width is much larger. Thus, in the small K 
regime, we have a much longer range interaction than in the large K regime. 
This can also be seen from the averaged correlation function C{\i—j\) = {Cij), 
where (• • •) denotes the ensemble average over particles. We have investigated 
this function for a wide range of K. We tried to use different forms to best- 
fit this function, but found that it is very hard to see whether the decay is 
exponential, power- law or logarithmic. It seems that it obeys different laws in 
different ranges. Nevertheless, the correlation function in the large K regime 
decays faster than that in the small K regime, as illustrated in Fig. 10 for 
K — 0.5 and K — b, respectively. 

On the other hand, due to discreteness, it is also very difficult to define a 
parameter such as the width of the correlation function. As an alternative, 
we chose C(l), which describes the correlation between the nearest-neighbor 
particles. In Fig. 11, we plot C~^(l) against K on a logarithmic scale. The 
signature of the transition is very similar to Fig. 8. In the small K regime, 
C(l) is a constant (about 0.5); and in the large K regime, it decreases with 
K. Since the transition is largely smeared by quantum fluctuations, one can't 
define a critical Kc as in the classical case. 



Cii 




(31) 
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We know that if the correlation function decays exponentially in the short 
range, then InC^^(l) = /^^, which is just the Lyapunov exponent a in the 
classical case [2], / being the correlation length. Thus, the parameter C(l) 
captures the short-range correlation. Fig. 11 mimics the classical plot of the 
correlation length proposed by Aubry [2]. In the classical case, for K < Kc, 
I — oo and (7 = 0, whereas for K > Kc, I is finite and a > 0. However, in the 
quantum FK model, we do not have infinitely long interaction in the small K 
regime. Instead, we have a finite interaction range, but it is much longer than 
that in the large K regime. 



5 Concluding remcirks 

The Frenkel-Kontorova model is a simple model. However, despite its decep- 
tively simple appearance, it exhibits very complex behavior and contains very 
rich physics. It has found applications in a wide range of fields, from nonlin- 
ear dynamics to condensed matter physics to tribology. However, most of the 
studies have so far been confined to the classical Frenkel-Kontorova model. 
Very little work has been done on the quantum Frenkel-Kontorova model. We 
think this is indeed quite regrettable. The quantum Frenkel-Kontorova model 
is not only interesting from a theoretical point of view but is also important 
from a practical point of view. The work reported in this review represents 
only a first step in this direction. We however hope that it will serve as an im- 
petus to further studies. It is our conviction that a thorough understanding of 
the quantum Frenkel-Kontorova model will lead to important breakthroughs 
both in theory and applications. 
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hl-figl 




Fig. 1. The classical hull function and ^-function in the subcritical regime {K = 0.5) 
and the supercritical regime {K = 5). In (b) and (d), the solid line is sinX^. 



U-fig2 




Fig. 2. The quantum hull function and ^-function from QMC, h = 0.2. The solid 
line in (b) and (d) is sinX. 
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hl-fig3 




Fig. 3. AXf = {{Xi — XiY) = hGii, calculated from QMC. The parameters are the 
same as in Fig. 2. 



M-fig4 




Fig. 4. hGij as a function of i and j. The data are obtained from QMC with h = 0.2, 
K = 5. 
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Fig. 9. Contour plot of the correlation function Cij for K = 0.5 and K = 5 at 
h = l. 
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Fig. 10. {Cij) as a function of \i — j\. The average is taken over the range of i and 

j, n = i. 
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Fig. 11. C ^(1) as a function of if for ^ = 1 and 2, respectively. 
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